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Oh' Abstract. 

The competition between the drive and stabilization of plasma microinstabilities 
by sheared flow is investigated, focusing on the ion temperature gradient mode, 
r^ Using a twisting mode representation in sheared slab geometry, the characteristic 

Cu , equations have been formulated for a dissipative fluid model, developed rigorously 

from the gyrokinetic equation. They clearly show that perpendicular flow 
shear convects perturbations along the field at a speed we denote by Mc s 
^fl ' (where c s is the sound speed), whilst parallel flow shear enters as an instability 

O ' driving term analogous to the usual temperature and density gradient effects. 

' Ty!. , For sufficiently strong perpendicular flow shear, M > 1, the propagation of 

^Li the system characteristics is unidirectional and no unstable eigenmodes may 

(— I form. Perturbations are swept along the field, to be ultimately dissipated as 

they are sheared ever more strongly. Numerical studies of the equations also 

reveal the existence of stable regions when M < 1, where the driving terms 

conflict. However, in both cases transitory perturbations exist, which could attain 

f%l , substantial amplitudes before decaying. Indeed, for M 2> 1, they are shown to 

^. i exponentiate vM times. This may provide a subcritical route to turbulence in 

f""s ■ tokamaks. 
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PACS numbers: 52.30.Gz, 52.35. Qz, 52.35.Ra 
1. Introduction 



The confinement of heat and particles in tokamaks is determined by fine scale 
turbulence. This turbulence is driven by the equilibrium gradients in temperature, 
^^ , density and flow velocity pQ. A central goal of fusion research is to understand and 

suppress this turbulence. There is increasing evidence that tokamak turbulence is 
strongly affected by sheared flows. For example, recent experiments on DIII-D show 
confinement increasing with toroidal flow shear [2]. Furthermore transport barriers, 
thin regions of the plasma with very large temperature gradients, have strongly sheared 
flow velocities j3j. It is presumed, but perhaps not proven, that turbulence and 
therefore heat diffusivity is suppressed by the sheared flows in the transport barriers. 
There is considerable literature on the stabilization of plasma instabilities by sheared 
flows (sec for example [IHE])- More generally it has been argued [9] that flow shear 
suppresses turbulence by shearing apart the turbulent eddies. There is also intense 
interest in the regulation of plasma turbulence by self-generated sheared flows — see 
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for example |10H13] - - and the transport of plasma angular momentum, which is 
required in order to understand self-consistent flow profiles — see for example |14H16] . 

Of course, not all shear flows are stabilizing. Indeed shear flows parallel to the 
magnetic field are destabilizing (this has been known since the 1960s, see [1711X5] ). 
Also perpendicular shear flows can lower the instability threshold of gravitationally 
driven MHD instabilities — see [19] . Typically strong flows in tokamaks are toroidal 
and are thus a mixture of parallel and perpendicular flow. What then is the optimum 
shear flow for confinement? Is increasing toroidal flow shear always desirable? In this 
paper we begin an examination of these questions by investigating the effect of flow 
shear on the ion temperature instability in a fluid slab model. Although more realistic 
descriptions have been investigated numerically [7ll8l l20H22] . this simple model allows 
a clear demonstration of some important features. 

The dynamics of the instability is investigated in coordinates that simultaneously 
shear with the flow in time and with the magnetic field in space. These transformations 
individually have a long history (see [231E3]) and have been used in combination more 
recently [19,25 27. Instabilities then take the form of twisting eddies that travel 
along the field at a speed Uf (sec figure Q] and equation ((7j for a definition). We define 
a Mach number M which is the ratio of u/ to the sound speed. Whilst others have 
used the more familiar Fourier representation [6l (HI EH1 E9j , the "twisting-shearing" 
representation used here clarifies aspects of the physics. The fluid approximation 
analyzed here does not, of course, treat collisionless dissipation (Landau damping), 
nor does it accurately capture ion finite Larmor radius (FLR) effects and therefore 
cannot model accurately modern fusion devices. However it does provide a simple 
model to investigate the physics of sheared flow. 

We focus on the linear evolution of the ion temperature gradient mode (referred 
to as the ITG) and the mode driven by the parallel velocity gradient /shear (the 
PVG). Note that, as has been discussed previously [6j[28l[29], the PVG mode differs 
from the typical incompressible Kelvin Hclmholtz instability, the free energy source 
being shear in the flow parallel to the magnetic field and compressibility effects being 
essential, allowing destabilization of sound waves. We show that the twisting mode is 
localized where the perpendicular gradients are small so as to minimize the collisional 
dissipation. The exponential stability for moderate perpendicular flow shear (M < 1) 
is investigated numerically. 

At large sheared flow, M > 1, exponential growth is stabilized as the 
perturbations are swept downstream. This result is proved in section [5] by considering 
the system characteristics. In an important paper, Waelbroeck et al. [5] showed 
(using the Fourier representation) that the parallel velocity gradient mode in a cold 
ion plasma is stabilized above a critical perpendicular flow shear. Their result can 
be translated into the criterion M > 1 for their equations in the twisting shearing 
representation. A similar result was observed by Howes et al. [19], who showed that 
the MHD gravity-driven instability of a sheared slab is stabilized when u/ exceeds the 
Alfven speed. 

It is well known that sustained turbulence is possible in systems without linear 
instability |30H32| . Indeed it is a common cause of turbulence in fluid systems such as 
plane channel flow [33] . Waelbroeck et al. [34) suggested that transient perturbations 
driven by parallel flow shear could lead to turbulence above the critical perpendicular 
flow shear. We present in section 15.21 an asymptotic calculation of transient growth 
in the large shear flow limit M 3> 1. This solution demonstrates that for M ^> 1 
transient perturbations exponentiate many times (\^M times) before decaying. We 
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have not considered the non-linear implications of these transient perturbations — but 
one expects that they amplify enough to sustain turbulence. In fact recent toroidal 
gyro-kinetic simulations with sheared flow by Barnes et al. [3Tj and Highcock et al. [2"2"] 
have identified regimes where exponential instabilities are stable but transient growth 
triggers turbulence. 

The analysis is presented as follows. In section [2j the collisional plasma response 
to an electrostatic perturbation in a magnetized slab is summarized, highlighting the 
distinct effects of parallel and perpendicular flow shear. The linear evolution of the 
system is considered for the remainder of the paper. The local dispersion relation 
(ignoring collisional dissipation) is discussed briefly in section [3j focussing on the 
ITG and PVG modes. Then the effect of collisional dissipation and the localization 
by a sheared field without flow is considered in section SJ In section [5J the evolution 
equations of the system are cast into characteristic form and the behaviour of transient 
perturbations analysed in the large shear flow limit, M ^> 1. Results from a numerical 
study of the system are presented in section [51 supporting the analytical results 
developed and extending the investigation, particularly to moderate values of M. 
Finally, a discussion of the implications of our results and a summary of the work are 
given in the conclusions, section [71 

2. System Equations 

In this section we present the system of equations that we use to analyze the stability 
of ion temperature gradient and velocity gradient driven modes. 

2.1. Geometry 

In order to highlight the main physical effects of sheared flow, we consider the 
stability of a simplified geometry — the infinite plasma slab in the presence of sheared 
background magnetic and flow fields of the form: 

B = B ° ( 2 + T/) ' (1) 

V = Vb^-e v , (2) 

where e v is a unit vector in the (y, z)-plane and l s , L v are, respectively, the 
characteristic scale length of the magnetic field and flow variation in the x-direction. 
This geometry is illustrated in figure [IJ which is discussed further below. The variation 
of the background plasma density and temperature, which drives instability, is also 
taken to be in the x-direction (see section |2~2|) . This system has been analyzed many 
times, including studies of the stability of the ion temperature gradient driven mode 
with sheared background flow [7j[28]. Note that we have chosen a frame in which 
there is no flow on the x — surface. The instabilities are assumed to be localized 
where (x/L v , x/l s ) <C 1 so (TTJ) and © can be thought of as Taylor expansions about 
x = of more general field and flow profiles. The sheared flow has components parallel 
and perpendicular to the magnetic field. As we show, these two components produce 
distinctly different physical effects. 

To analyze the system, it is most useful to implement a double-shearing coordinate 
transformation, which was outlined for a slab in |19j . Such a transformation was 
introduced into the torus by Cooper [35] (see also [20,25,26). This removes the 
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complications of both the magnetic field shear and the time dependence associated 
with the equilibrium flow. It also aligns coordinate lines with the characteristics of the 
plasma response (the guiding centre motion and sound waves) and therefore simplifies 
both the analysis and the interpretation. In order to treat the effect of magnetic 
field shear, we first make the transformation [24] to the set of field line following 
coordinates: 

(3) 
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The component of the equilibrium flow perpendicular to the magnetic field 
convects (see figure |Tfb)) the magnetic field structure, the guiding centres and the 
sound waves. The time dependent transformation 

y*=y-V a ^-(e v -y)t (5) 

follows the shearing arising from the y (perpendicular) flow component. This type 
of transformation was first introduced by Kelvin |23| to analyze instabilities due to 
sheared flow in fluids. Combining the two transformations yields the final coordinate 
transformation |19j : 

t = t, z = z + Uft, y = y — —z , x = x, (6) 



1, 



where we define the velocity: 



Uf = V Q -^e v -y. (7) 

Thus we have transformed to a moving frame travelling along the z direction 
(essentially the field lines) at velocity — Uf. In this frame the combined transformation 
is identical to © and therefore independent of time. The coordinates yield the 
simplifications: 

B ■ Vx 1 = B ■ Vy' = 0, B V = B -^, (8) 

oz 

d , r ~ d d d ,„. 

di + V °- y dy- = dF +Uf d?- (9) 

Note that in the moving frame the plasma flows along the field line at velocity Uf. 
Perturbations in this system take the form of eddies which are localized in x and 
extended along the field line — cross sections through a single eddy (red ovals) are 
shown at various z positions in figure []Ja). The eddy tilts to remain on a surface 
of constant y' , along which the plasma has its characteristic response. With time 
this surface, and hence the eddy, is twisted by the perpendicular component of the 
sheared flow, sec ©. This is shown in figure HJb) and the drive aligned position of the 
eddy (eddy lies parallel to x) retreats along the field at the speed Uf. For sufficiently 
large Uf the eddy dynamics is too slow to allow this response and the mode structure 
indicated becomes transitory (see section [5]). 
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Figure 1: Geometry of the shearing magnetic field B and background flow Vq. Looking 
in the — z direction, the field is represented by flux tubes, twisting above and below 
the plane of x = 0. Dashed grids are a guide, both l s and L v are taken to be negative. 
(a) Cross sections (red ovals) through a typical eddy at an initial time t = 0. (b) 
Flux tubes are advected along y by the perpendicular component of the sheared flow, 
see ©. The eddy is thus twisted and its drive aligned position (eddy parallel to x) 
retreats along the field at speed Uf = V (l s /L v ) e v • y. 



2.2. Plasma Response 

The equations used to describe the plasma response will now be presented. The 
plasma is taken to be quasineutral, consisting of electrons and a single (hydrogenic) ion 
species, of charge e and mass m. The equilibrium density of both species is hq. Only 
electrostatic perturbations are considered, with the perturbation frequency, u, much 
less than the ion cyclotron frequency, £1 = eBo/m. For simplicity, the equilibrium ion 
and electron temperatures are taken to be equal in the final form of the equations and 
denoted by T. We consider instabilities with phase velocities comparable to the ion 
thermal speed and therefore much slower than the electrons. The electrons thus have 
time to set up thermal equilibrium and so are taken to have an isothermal Boltzmann 
response to the perturbing potential <f>, giving the perturbed electron density 

(10) 



6n e _ = 



T 



n . 



<1, 



(11) 



The ions are described by collisional fluid equations; the full details of their 
derivation are given in the appendix. The equations are derived from the gyro-kinetic 
equation [SHIES] thus they embody the usual gyro-kinetic ordering: 

n \k±J \f Q/ 

as a primary expansion. Parallel (or perpendicular) is taken with respect to the 
equilibrium magnetic field, k is the wavenumbcr of the perturbation and Sf is the 
perturbed distribution function, where Fq represents the bulk distribution. We make 
a subsidiary collisional expansion, as shown in the appendix, with: 

v 3> u), oj* , w t hfc||, Ufk\\, vk p . (12) 

Here v is the ion-ion collision frequency, uj* represents the drift frequencies associated 
with the background gradients (see section [231) . v t h = (2T/m) is the ion thermal 
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velocity and p = v^/fl, is the ion gyroradius. All quantities on the right hand side of 
(|12[) are treated as the same order in the derivation of the fluid equations. Note that 

kp~QLjujv\ < 1. 

Upon expanding the gyro-kinetic equation with respect to uj/v, wc obtain to 
lowest order a perturbed Maxwellian distribution (see appendix). The perturbed 
density, Sn, parallel velocity, SVu , and temperature, ST, of this Maxwellian obey 
fluid equations — the moment equations. The three moment equations for particle, 
parallel momentum and energy conservation may be combined with the requirement 
of quasineutrality on the perturbed electron and ion densities. This gives a 
closed set of equations for the evolution of the perturbed variables. We employ 
the following normalizations, which are used in many gyro-kinetic codes (see for 
example 38 ). They indicate the typical features associated with the unstable 
modes: the rapid perpendicular and long parallel spatial dependence, the characteristic 
acoustic timescale and the amplitude scaling. 

p s y, z' = l s z, t = —i, (13) 

Cs 

ST l s 8nl s 

^r — , n= , (14) 

i Ps n p s 

where p s = c s /Q is the sound Larmor radius associated with the sound speed 
c s = y/i^e + li) T/m. In this collisional model, j e = 1 (isothermal) and -fi = 5/3 
(adiabatic). The background gradients provide the instability drive. The scale- lengths 
for the equilibrium density, parallel velocity, temperature and the Mach number 
associated with the moving frame are: 

Id, Id, 

l n dx It dx 

1 1 Vn u f 

T = - -e v -z, M=-i-. (15) 

i>v ^v C s Cs 

The tildes denoting the final transformations will henceforth be dropped for 
convenience. 

The set of three equations describing the evolution of a perturbation of the system 
in the doubly sheared system of coordinates thus takes the form: 

d , r d \ dV 3 ; s dn 

oi + M oi) n+ ^ = -^W <16) 

^"hir + m^ + lM-lZfrrtv, (17) 

where 

V 2 — ( — — 
dy 2 \ dx dy 

and with the unit vector defined as b = B/B, for two arbitrary functions p and q the 
normalized Poisson bracket form is: 

{M}=(VpxV ? ).b4^^. (20) 

ox ay ay ox 



z ~ ) (19) 
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The time derivative appearing in these equations has a convective form: -Sr + M ■§- 
- in the moving frame the plasma is flowing along the field lines at a normalized 
velocity M. As discussed above, the potential <j) is proportional to n (see ((TU|)), so 
the perturbed E x B drift in our normalized time and spatial coordinates is given 
by Vn x b. Thus, nonlinear terms describe convection of the perturbations by the 
perturbed E x B drift. Since {n, n} = 0, there is no nonlincarity in the density 
evolution equation. The instability drive comes from the first term on the right 
hand side of each of equations (fTMTSj) . They are, respectively, the convection of 
the equilibrium density, parallel velocity and temperature along their gradients by the 
perturbed Ex B drift. Dissipation results from the diffusive viscosity and conductivity 
terms. These are derived in the appendix and have the following normalized forms: 

/9 1\ f2 l s n e 4 \nA 
giving a Prandtl number of: 



40' 4/ V 3 87T 3 / 2 e 2 T 2 



Pr = ^ = 4 (22) 

X± 10 
In the sheared coordinates the dissipative term becomes large at large z — this is 
just the mathematical manifestation of the shearing of the (field) coordinate lines. 
However, as we shall see, the perturbations follow the (field) coordinate lines and 
therefore the effective dissipation does indeed increase along field lines. As ion-ion 
collisions conserve ion momentum and therefore do not produce particle transport, 
there is no particle diffusion in the density equation. If we had included ion-electron 
collisions we would have obtained a particle diffusion of order: D± ~ yjm e /mi \±- 
Note also the limitations pointed out in the introduction: a fluid model such as this 
cannot capture collisionless dissipation or FLR effects. However, its particular strength 
lies in the clear representation of the physical effects of sheared flow which it allows. 
The background flow shear enters the system of equations ([T6l - fT8|) in two forms. 
First, as a parallel (along z) convective velocity, M, due to the perpendicular 
component of the flow shear, which acts on all the perturbed variables. Second as 
an instability drive, oc l s /l v , from the parallel component of the flow shear, which 
only acts on the parallel velocity perturbation. The magnitudes of these flow effects 
arc linked, as both are proportional to the magnitude, Vo, of the imposed background 
flow, so: 

M e v • y . , 

t, Tj-r = r=tan^„. (23) 

yl>s/ tv) *^v ' z 

The aim of this paper is to understand the interaction of these two effects and whether 
sufficient flow shear can stabilize the system via strong convection, which will be seen 
to sweep the perturbations along the field lines into the dissipative regions at large z. 

2.3. Linearized System 

For the remainder of this paper we neglect the nonlinear terms in the equation set 
above and restrict ourselves to the investigation of the linear evolution of the system. 
All fields arc taken to vary as cxp(iky) multiplied by a function of z. The perpendicular 
gradient operator then becomes: V^ = — k 2 (f + z 2 ) and we define: v^ = k 2 v± and 
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Xk = k 2 x±- The linearized set of equations then takes the form: 



9 M d> 

dt dz / 


jn+ — = iuj n n, 




(24) 


d Jr d\ 

d-t +M dzj 


)V + l§- z (2n + T): 


= iu>*n - v k (1 + z 2 ) V, 


(25) 


d-t +M dzj 


1T 23V . , 


- Xfe (1 + ^ 2 ) T. 


(26) 



The following effective drift frequencies have been introduced, characterizing the three 
background driving gradients: 

ok l s ok l s ok l s /r>~7\ 

o L n o i v o It 

For a perturbation of the form cxp(iky + ik x x) the equations are identical if we make 
the transformation z —} (z — k x /k). As noted above the dissipative terms increase 
strongly with distance along the field line. 

In the following two sections analytical limits of this third order system arc 
considered, which highlight the effects of parallel and perpendicular flow shear on 
the system. We consider the local dispersion relation first. 

3. Local Dispersion Relation 

The dynamics of the basic modes present in this third-order system may be identified 
by considering the local, dissipationless form of the above equations. By neglecting the 
effects of dissipation, v k =0, Xfc — 0i the terms explicitly proportional to z vanish, so 
we may look for solutions in the form exp(— iut + ikuz). Noting the purely convective 
effect of the perpendicular flow shear, M , we define the frequency in the laboratory 
frame as: 

u' =u-k\\M. (28) 

The dispersion relation is then: 

lo' 2 ( W ' + <) - W '*|| (fey - u>*) + -1 (|4 - w;J = 0. (29) 

When w* = oj* = we recover the usual cubic dispersion relation for the 
collisional ion temperature gradient (ITG) instability |39j . The same limit for 
instability, cj t > (16/9V3) k\\, applies here with finite M, but the corresponding- 
oscillatory frequency of the mode has changed, by — k\\M. For ui T ^> feii, ((2^)) reduces 
to the original ITG dispersion relation: u>' 3 = — 3k 2 u) T /8, derived by Rudakov and 
Sagdeev gD]. 

In the limit w* = lo t = the dispersion relation simplifies to a quadratic: 

J 2 = fc|| (fc|| - u*) . (30) 

This describes the parallel velocity gradient (PVG) instability introduced by previous 
authors [BlflS]. with instability if: 

<>k\\- ( 31 ) 

The most unstable mode has k 1 J iax = uj*/2 and thus a growth rate r y m ax = w*/2. The 
stability of the o;*-driven mode is seen from (|31[) to be sensitive to the relative direction 
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of the mode propagation and the flow shear, specifically via the ratio kj (kul v ). A 
clear discussion of the dynamics of this parallel flow shear-driven mode was given by 
Catto et al. [IB], among others. Perturbing the density causes a perturbed electric 
field, via Sn e = (e(f>/T)no. The E x B motion in the perturbed electric field, with 
direction defined by fc, convects the background parallel flow. The perturbed flow 
caused by this convection enhances (for k/k\\l v > 0) or opposes (for k/k^l v < 0) the 
parallel motion causing the density perturbation. Instability results when the rate of 
increase of the flow by the convection overcomes the deceleration due to the parallel 
pressure gradient — this happens when u>* > fc|| . 

The parallel and perpendicular flow shear are related, as noted in section [2721 so 
the stability criterion pip may be translated into a requirement on the Mach number. 
We define: 

wl = aM, a = -A^-4, (32) 

8 e v • y 

and the angles 9 and 9 V , giving the direction of mode propagation and of the 
background flow, with respect to the magnetic field: 

tan 9 = — , tan 9 V = — — - . (33) 

fcii e v • z 

Thus (f3~Tj) implies that, for instability: 

aM>k\\. (34) 

In the case when both uj* and M are positive, we obtain the instructive rearrangement: 

M>M 8tan^ 
a A tant> 
(Note that we may always take M > without loss of generality due to the symmetry 
of the system.) This form will prove convenient in section [5] and simply indicates that 
for a background flow velocity with a finite parallel component, modes with a finite 
parallel wavevector require a stronger drive to become unstable. Note that the mode 
frequency is given by (J3T)|) . 

Various stability boundaries may be derived from (|29p. It is illuminating to 
consider the two-dimensional stability boundary in the (u>*,u>^) -plane, which shows 
the impact of parallel velocity shear on the ITG instability. The dispersion relation 
in this case reduces to: 

uj' 3 - u/fc|| (fc|| - uj*) + -fcjju4 = 0. (36) 



Instability exists if either: 
u 

J- 



" > 1 (37) 



or 



16 



3/2 



>^-(l-^f\ . (38) 

This can be seen to reduce to the ITG stability criterion for o>* = 0. The stability 
boundary is shown in figure [2] and the asymmetric effect of background parallel flow 
shear on stability may be clearly seen. 
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Figure 2: Interaction of the ITG and PVG modes in the local dissipationless limit. 
Contours of the growth rate (divided by k\\) of the fastest growing solution to (|36|) arc 
shown, with u>* = 0. 



4. Parallel Flow Instability with M = 

The dissipativc terms which were neglected in the previous section will act to 
damp the perturbations at large distances along the field line. Analytical solutions 
demonstrating this may be obtained by considering the case with no perpendicular flow 
shear, AI = 0. For tractability we retain only viscous dissipation in the system. (As 
may be seen from equations (|24H26|) . this corresponds to a system with an artificially 
large Prandtl number, vxjxx ^ 1-) The linearized equations then reduce to: 

dn dV 

■wr + -z- = iu n n, 

at az 

' (2n + T) = iu* v n - v h (l + z 2 ) V, (40) 



(39) 



dV 
~dt 


3 d 

+ 8dz' 


dT 
~dt 


2dV 
+ 3~dz~ 



(41) 



No explicit time dependence appears, so we may look for a solution which varies as 
exp(7<). The set can then be formed into a single mode equation, written here for 
the perturbed parallel velocity: 



d bj* \ d a , on 

«- " ? '-r )1T V ~ I (7 + "k + v k z 2 ) v = o. 
az b az b v ' 



(42) 



The effect of the drives wj and w* are given by the parameters a and b, defined as 
a = 7 - iw„, 

b = H 

4 7 \2 

It is convenient to define: 
V 



(43) 
(44) 



s (z) exp (iuj*z/2b) exp (jt) , 

Note that s (z) must converge sufficiently rapidly as 
\z\ — > 00. Equation P2^1 reduces to the form: 



dz 2 



<V a a 2 

-w-~b {l+Vk) -r kZ 



(45) 
— > 00 to ensure that V — > as 

: 0. (46) 
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Bound state solutions (those satisfying the boundary conditions V — >• as \z\ — >• oo) 
require that the real part of y/a/b be positive. Upon defining the variable: 

the mode equation may be cast into quantum harmonic oscillator form: 
Fs 

The parallel structure of the mode, described by s (z), therefore has the form: 



, [e-e]s = 0- (48) 



s (z) = H„ 



-yjav k /b (z 2 /2) 



(49) 



where H n are the Hermite polynomials. The mode thus corresponds to waves, of n 
nodes in the parallel direction, trapped between dissipative regions far along the field 
line. There is also the overlying rapid oscillation, given by (|45l) . so: 

1/4 



V = H n 



ffi 



i/XT + '^l 



The quantized eigenvalues in (|48l) . e = In + 1, for positive integer n > 0, 
correspond to the allowed growth rates, such that: 

£( 7 + I *) + (2„+l)^-^=0. (51) 

In the limit w^ = w* = 0, the effect of dissipation on the parallel flow shear driven 
mode is seen to produce a growth rate satisfying: 

,*2 

,2 



7 2 + ^7+(2n+l)>ia=-^. (52) 

Thus a growing mode exists for all finite values of u>* and the growth rate is increased 
for a given drive strength if the dissipation, i>k, is reduced. When V). = 0, we recover 
the relation 7 = w*/2 of section|3l At large k (y^ 3> w*) the mode is still growing but 
dominated by viscosity 7 ~ uj* 2 /4vk- In reality large k modes will be damped cither 
by particle transport or finite Larmor radius corrections — effects excluded from this 
analysis by our ordering. 

Shown in figure [3] is a validation of equations (|51[) and (|45p against a numerical 
solution of the full set of equations (|2~4T[26| in the limit M = 0, \ = 0. (A description 
of the code is given in section [SJ where results of numerical investigations of the 
full linearized system defined in section 12.31 are presented.) The numerical solution 
displayed is the fastest growing mode found, whilst the analytical solution plotted 
corresponds to the case n = 0. 

5. Impact of Convection 

As was seen in section [521 the perpendicular flow shear acts to convect perturbations 
along the system. Dissipation, aided by the strong magnetic field shear, acts to remove 
all but the density perturbation at large distances along the field lines. Thus we may 
anticipate that rapid convection will sweep growing instabilities into the dissipative 
region and force their decay. This is motivated in the following section by considering 
the characteristics of the linear system, equations (|24fj2"oT) . In such cases, the system 
would be considered linearly stable, as no exponentially growing cigenmodes would be 
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Figure 3: (a) Growth rate as a function of k. Squares are obtained from the numerical 
solution of the system of equations (J24l - |2"6"|) . with parameters: l s /l n = 3, I s /It = 6, 
l s /lv = 9, M = 0, v = 3 and x — 0- The corresponding analytical solution (|5"Tj) for 
n = is also shown (solid line), (b) The numerically obtained velocity eigenfunction 
(absolute value) for k = 0.3 (dotted line) with the analytical expression (|50]) (solid 
line). 



found. However, an unstable mode would be able to grow for a finite time. The form 
of such transitory solutions, which begin to grow exponentially but are then forced to 
decay, is determined in section IS31 As discussed in the Introduction, it is of interest 
to determine if transitory solutions grow to sufficient amplitude to trigger non-linear 
effects. 



5.1. Characteristic Equations 

The linear set of equations describing the system, equations 
conveniently written in characteristics form. Defining: 





s 


= -T-n, C ± = V±-(n+- 
2 4 V 2 


we obtain: 


ot oz 


s = ^(c+- C --§)-±» (! + **) 


! + (M ± 
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c± =4( 0+ -c--l) 
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\u k (1 + z 2 ) (C+ + C-) T \ Xk (1 + z 2 ) (V -C~ + ^- 



can be 

(53) 

(54) 
(55) 



where 



(56) 

These characteristics represent three waves: the entropy mode propagating the 
perturbed specific entropy, S, at speed M, the forward propagating sound wave 
propagating C + at speed M+l and the backward propagating sound wave propagating 
C~ at speed M — 1. The terms on the right hand sides of ((55J) couple the three waves 
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but do not change their propagation speeds. An initial perturbation which is localized 
in z between z = a and z = b at time t = (i.e. the function describing the initial 
perturbation has compact support in z) must be localized between z = a + {AI — l)£i 
and z — b + (AI + l)ti at time t = t\. Eigcnmodes can form when M < 1, by the 
combination of oppositely travelling waves coupled by the right hand sides of (j55j) . 

As M is increased, we see from (|55p that the speed of the forward travelling 
characteristic is enhanced, whilst that of the backward is reduced. Thus, when M > 1 
is reached, all of the characteristics of the system are forward going, due to the 
convective effect of the perpendicular shear of the background flow. In this case 
no eigenmode may be formed in the system. Indeed the initial perturbation is swept 
forward, since z = a + (AI — l)ti and z — b + (AI + l)ti both increase with the time 
t\. This discussion is borne out by the numerical solutions presented in section O 
see figure [5J As described in section 21 velocity and temperature perturbations will 
damp at large distances along the field line. Note that the density perturbation will 
remain, due to the form of (|!M]) . as will also be seen in section[5J Therefore, for M > 1, 
any initial unstable perturbation will be swept into the dissipative region and forced 
to decay. Whilst such a perturbation can grow exponentially for a finite time, its 
behaviour would not be captured by a traditional eigenvalue analysis. In the following 
section we consider the evolution of such a transitory perturbation. 

5.2. Transitory Solution 

In order to determine the behaviour of transitory instabilities, we consider a system 
with strong shear of the perpendicular background flow, that is M > 1. An analytic 
solution may be obtained upon retaining only the driving term due to the parallel flow 
shear, which is associated with the strong convection via (f3"2j) . and formally neglecting 
the effect of thermal conductivity, as in section 0] The solution in this limit captures 
the nature of the transitory behaviour, which arises due to the strong convection 
rather than being dependent on the details of the dissipation mechanism. Therefore we 
consider in this section the transitory evolution of the w*-driven instability described 
in section [3J The transitory behaviour of the full system was investigated numerically 
and the results are presented in section [5J 

It is convenient to analyse the system in the moving frame, using the coordinates 
(zo, t), where zq = z — Alt. With wj = a;* = Xk = 0, the system of equations (|24][2G| 
may be reduced to a single differential equation, written here in terms of the parallel 
velocity perturbation: 

d 2 V 8 2 V . ,dV d r, 2W1 

^ = ^-^-^ 1 + ^- (57) 

The time derivative is to be taken at constant zq. The rapid parallel oscillation 
inherent in this instability, which was present when Ad — (see section 2]) , may be 
removed upon substituting: V (z , t) = V (zq, t) exp [iuj*Zo/2}. This reduces ([57} to: 



d 2 V d 2 V uf~ 
dt 2 8zq 4 



dV 

[1 + zl + 2z Mt + M 2 t 2 ) — + 2AI {z + Mt) V 



(58) 



We anticipate that the initial exponential growth of the instability in time, 
associated here with the explicit w* term as seen in section [3j will be overcome after 
a finite time by dissipative effects when M S> 1. This exponential behaviour may 
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be conveniently characterized by introducing the time dependent "growth rate" 7 (i) , 
such that: V (zo,t) — V (zo,t) exp J j(t')dt' , satisfying: 



d 2 V 



dV 



W + 2l 'dt 



+ 7 V = 



d 2 V 
dzl 



-v - ^ 



(1 + z 2 + 2z Mt + M 2 t 2 ) Uv + ^ J + 2M (z + Mi) V 



(59) 



This evolution equation may be expanded in the large parameter, M. We seek a 
solution where: 



dV 



and ^ ~ O ( 



(60) 



Remembering that for short times 7 (£) will approximate the local solution of section[31 
where 7 ~ uj* oc M , the leading terms of ([5T))) arc of O (M 2 ) and yield: 
,*2 



a.' 



7 



=J- - ^M 2 t 2 7 . 



Hence the function 7 is given by: 
To 



= -T 2 + VV 



1, 




Mi, 



(61) 



(62) 



where we have taken the positive root of the radical, as the behaviour of the fastest 
growing instability is of most interest, and we have normalized to the local value of 
the growth rate in the case of instability, 70 = oj*/2 (see section [3]). For long times, 
7/70 —* (2t 2 ) -1 and exponential growth will no longer dominate the evolution of the 
perturbation. Therefore we may identify that the instability will grow exponentially 
for a characteristic time t c = -^70/ (^/UkM\ oc 1/vM, during which time the initial 

where a 



perturbation will be amplified by a factor of order exp [7oi c ] ~ exp ayM 

is independent of M. Clearly at large velocity shear (M) the transitory perturbations 
can grow significantly. 

The amplitude of the perturbation will turnover around t c and begin to decay 
due to dissipation. Thus the function V will change on the timescale r c (slower than 
the growth time 7^ 1 ) as anticipated in ([60)) . The behaviour of V will be governed by 
the next significant terms in ([59]) . which are of O (M 3 ' 2 ): 



2" 



dV 
~dt 



■v k Mt 



Alt 



dV 
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+ 2z -yV + 2MV 



(63) 



This is simply a first order equation for the time dependence of V, in which zq is a 
parameter, therefore with 7 and r given by ([62]): 



V(z ,t) = V(z ,i)exp 



V(z ,t) 



'1 



£irV 
M , 



vT 



exp 



d(r 2 ) 
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(64) 



,-1? 



where we have used the substitution r 2 = sinh -d and noted that sinh -d — cosh i3 = 
to obtain the final form. The function V has been introduced to capture the remaining 
finite spatial and temporal dependence of the initial perturbation amplitude, present 
in (El- 
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The perturbed parallel velocity for M ^ 1 thus takes the final form: 



V (z , t) = V (z , t) (r 2 - Vl + r 4 ) 



exp 



^o + ^7(<) + 



(65) 



In the initial stage of the instability, as M 3> 1 the final term in the exponential will 
dominate the evolution, giving exponential growth of the perturbation amplitude with 
7 (t) w 7q. However, at times longer than the saturation time: t> t c ~ 1/y/M, 7 (i) 
varies as i -2 and the asymptotic form of V (zo, i) becomes: 



$ 



* 2 / 1 \ 3/2 / 1/2 

i%, + f°JIL-f^ + V ( 1.236- 



(66) 



The exponential factor therefore saturates and the amplitude evolution in the decaying 
phase is determined by the algebraic prcfactor, which decays as t~ 2 . Thus, not only can 
transitory perturbations reach significant amplitude during their initial exponential 
growth, they linger in the system, decaying in amplitude only algebraically with 
time. The density perturbations do not decay after an initial amplification and 
only effects outside this analysis (either particle transport or finite Larmor radius 
effects) cause their eventual decay. Note that the total number of e-foldings is: 
Q.ly / kM/v± (cot V ) ■ As expected, large parallel flow shear gives large transitory 
growth. 

6. Numerical Results 

The results of numerical investigations of the linearized system defined in section 12.31 
will now be presented. We integrate equations (f5~JH55|) numerically with a second- 
order accurate upwind scheme, which naturally captures the forward or backward 
propagation of the characteristics. The equations are solved in a box of size —Zi n f < 
z < Zi n f, where Zmj varies between 50 and 200 depending on the parameters. The 
typical resolution is Az = 0.1. Convergence tests were performed to establish that 
the size of the domain and the resolution employed were appropriate; a benchmark of 
the code against an analytical solution has been shown in figure [3J The perpendicular 
viscosity (|21|) is fixed at u± = 3.0 for all cases presented: this relatively large value 
of v± ensures the consistency of our orderings. The analytically defined instability 
drives (27) depend on the wavenumber k as does the growth rate. Therefore, in this 
section we present the stability results as functions of the physical instability drives, 
the gradient scale length ratios (l s /ln,ls/lT, ls/lv)- In figures showing the maximum 
growth rate for a given set of physical parameters, the growth rate has been maximised 
over the wavenumber k, using the well known Brent's method |41j . Finally, for clarity 
we remove the density gradient drive throughout this section, that is l s /l n = in all 
cases. The results of tests with finite values of the density gradient did not alter the 
conclusions presented here. 

In figure 3] we illustrate the effect of the parallel velocity gradient on the ITG 
instability. A contour plot of the growth rate of the most unstable mode is shown 
as a function of the temperature, I s /It, and parallel velocity shear, l s /l v , drives, 
at fixed convective velocity M = 0.3. Note that l s /l v = c s M cot 8 V , so the rr-axis 
represents variation of both the angle and the magnitude of the flow. As we have 
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taken M < 1, cigenmodes may form in the system and the growth rate is well defined 
— see section HTT1 As seen, for a fixed value of l s /l v , the temperature gradient has to 
be raised above a certain critical threshold for the ITG to be unstable. Conversely, we 
see that the minimum value of the velocity gradient required for instability shifts to 
the right on the plot as the temperature gradient is raised. This figure thus illustrates 
the complex interplay between the two different drives and highlights the fact that 
increasing one drive may have a stabilising effect on the system if the other drive is 
kept fixed. Beyond M — 1 no cigenmodes arc formed in the system and no growth 
rate may be determined in the usual way, as discussed in section 15.11 




Figure 4: Contour plot of the maximum growth rate as a function of temperature and 
parallel velocity shear drive strengths. M = 0.3 and l s /l n = 0. Solid lines indicate 
zero growth rate. 



The more realistic case of fixed flow angle V is exemplified in figure [5J where 
we show contour plots of the maximum growth rate as a function of I s /It and l s /l v . 
We depict two cases: figure [5ja) has 9 V = 2°, so the flow is nearly parallel to the 
magnetic field, which is representative of a conventional large aspect ratio tokamak, 
whilst figure EJb) corresponds to the case 9 V = 45°, which is more pertinent to the 
case of a Spherical tokamak. These plots illustrate particularly the stability threshold 
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Figure 5: Contour plots of the maximum growth rate as a function of the temperature 
gradient drive and M for l s /l n = and (a) l s /l v = 30M, corresponding to 6 V = 2°; 
(b) l s /l v = M, corresponding to 8 V = 45°. For M > 1 no eigenvalues can be found. 
Solid lines indicate zero growth rate. 
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appearing at M = 1, which has been discussed in section [SJ They also show clearly 
the same non-trivial features seen in figure 0] we see that stable, or weakly unstable 
regions can also be found for M < 1 in certain ranges of the drive strengths. Consider, 
for example, fixed M and thus fixed l s /l v . We see that as the ITG drive I s /It is 
increased, the maximum growth rate in the system (which corresponds to a pure PVG 
mode at I s /It = 0) actually decreases at first until a region of stability is reached; this 
arises because the ITG and the PVG instabilities compete, rather than reinforce each 
other. As I s /It is increased even further, ^max eventually starts growing again and 
the system gradually transitions to a pure ITG mode. 

Notice also that as I s /It increases, the M value required to enter the stability 
band varies and a region of instability does not necessarily reappear as M is increased 
further, depending on the angle of the flow. This is consistent with the result obtained 
by Waelbroeck et al. [BJ, who showed numerically that the value of perpendicular 
velocity shear (essentially the value M used here) required for stability increased as 
Is/It increased, for l s /l v in the region of AM. Note that the marginal stability curves 
presented in [6] are consistent with M < 1, upon recognising that the sound speed is 
defined there for cold ions and that the stability curves are determined for T e ^ Ti > 0. 
As kinetic effects were retained, the precise definition of the sound speed used here 
for M is somewhat too simplified to reproduce the stability boundary of [6] exactly. 

The evolution of the perturbed fields in the two distinct types of stability region 
either side of M = 1 is illustrated in figure [6] Taking l s /l v = 30M, as in figure EJa), 
and choosing I s /It = 30 and k = 0.1, the left column of figure [6] shows the case 
M = 0.3 which lies in the stability band of figure [SJa) , whilst the right column has 
M = 1.2. Even though both cases are formally stable, since there is no exponentially 
growing eigenmode, there is an amplification of the initial perturbation in both cases. 
This amplification is more significant in the case of M = 1.2, particularly for the 
density. Notice in particular that whereas the temperature and velocity fluctuations 
are eventually damped, the density perturbation is seen to remain in the system, 
as anticipated in section 15.11 due to the lack of explicit dissipation in (|24jl . For 
contrast, the perturbed fields in the case of an unstable eigenmode from figure [5ja), 
with M = 0.3 and I s /It = 45, are shown in figure [7] Here k = 0.3, corresponding to 
the fastest growing mode. 

7. Discussion and Conclusions 

In this paper we have investigated the stability of ion temperature gradient and 
parallel velocity gradient driven modes in a sheared magnetic field with parallel and 
perpendicular sheared flows. The instability takes the form of twisting modes that 
travel along the field lines at a velocity Uf, which is proportional to the perpendicular 
flow shear, see (|7|). Modes can be excited by both the ion temperature gradient 
and parallel shear flow — in some cases the drives compete and instability can be 
suppressed and in others they combine to enhance instability (see figures 2] and [5]) . 
We define the Mach number of the moving perturbation as M = Uf/c s where c s is 
the sound speed. In section [5] we show that when M > 1, perturbations are swept 
downstream until they are damped away — thus only transitory growth can take 
place for M > 1. At very large flow shear the transitory growth is substantial: the 
number of exponentiations before decay is proportional to \/M, see section 15.21 and 
equation (l66l) . 

While it would be foolish to conclude too much about the behaviour of fusion 
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Figure 6: From top to bottom: time history plots of density, temperature and velocity 
(absolute values) for l s /l n = 0, I s /It = 30, l s /l v = 30M and k = 0.1. Left column 
shows case M = 0.3 cf. figure Eta); right column shows case M = 1.2. Each field is 
normalized to its maximum amplitude at t = so colors show the amplification factor 
of the initial perturbation. 



devices from this simple model, a number of trends stimulate some speculation. First, 
one expects that when Uf exceeds the characteristic propagation of the instability, 
exponential growth is suppressed. The characteristic propagation obviously depends 
on the instability. For example, electron instabilities have propagation speeds of order 
of the electron thermal velocity. It also seems likely that, for ion temperature gradient 
driven modes, turbulent suppression is maximized at values of M close to, or just 
above, one. Increasing the flow shear above this value may lead to the destabilization 
of transient instabilities and strong resulting turbulence. Barnes et al. [21] and 
Highcock et al. [22| have recently shown that this trend is indeed seen in gyro-kinetic 
simulations. Since turbulence is suppressed by perpendicular shear flow but enhanced 
by parallel flow we expect the sheared toroidal flow to give more suppression in regions 
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Figure 7: From left to right: time history plots of density, temperature and velocity 
(absolute values, log scale) for l s /l n — 0, I s /It = 45 and l s /l v = 30M cf. figure Ufa). 
Here k = k max = 0.3, which is the fastest growing mode (j ma x = 0.17). Only 
the central fraction of the simulation box shown for clarity (simulation domain is 
—200 < z < 200). Each field is normalized to its maximum amplitude at t = so 
colors show the amplification factor of the initial perturbation. 



of strong poloidal field. Indeed simulations by Roach et al. [42j confirm this trend. 
This also suggests strong suppression in Spherical tokamaks, where the poloidal field 
often exceeds the toroidal field on the outside of the magnetic surfaces (where the 
curvature drive is destabilizing). Work is in progress to establish more clearly the 
trends in shear flow stabilization and to optimize transport. 
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Appendix: Derivation of Moment Equations 

In this appendix we give details of the derivation of the set of moment equations 
presented in section 12.21 As discussed there, the electrons are taken to have an 
isothermal, or Boltzmann, density response, <5n e , to the perturbing potential, <fi: 

Sn e = n — , (67) 

where T e is the equilibrium electron temperature. The ion response is obtained from 
the gyro-kinetic equation [7ll36|l37] . The ion distribution function, /, correct to first 
order in the gyro-kinetic expansion in oj/Q, is: 

f = F (e t R) + h(R,w h wx,t), (68) 

Sf = -^p^F + h(R,w h wx,t), (69) 



Understanding the effect of sheared flow on microinstabilities 



20 



where 5f is the perturbation from equilibrium, the particle energy is e = mw 2 /2 + 
ecf>(r,t), the particle velocity is v, the velocity variable w = v — Vo and the 
guiding center position satisfies R = r — b x v/Vt = r — p. With the ion 
thermal velocity v 2 h = 2Ti/m for an equlibrium ion temperature Tj, the Maxwcllian 
F = (n /ir 3 / 2 vl) C xp(-e/TA. 

The distribution function of gyrocentcrs, h, is independent of the gyroangle of 
the particle motion, £, and is dchned by: 
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(70) 



The collision operator appearing here, C , is the linearized part of the self-collision 
operator [13] acting on the ion distribution function F + h and is discussed further 
below. The angled brackets denote the average of the enclosed quantity over the 
gyroangle at constant R: (A (r)) = (2tt)~ § A (R + p) dC, and the Poisson bracket is 
defined as follows, where the spatial gradient is taken at constant w: 



{(0),M = (v(0)xv/i)-b. 



(71) 



The ion response is driven by the background density and temperature gradients, as 
well as the parallel flow shear: 
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The ion sound speed, c s , is defined in section |2~2"1 The double shearing coordinate 
transformation discussed in section [2~T1 is now implemented, with Uf = Vq (l s /L v ) e v -y, 
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di = dt> + Uf Ih> ~ Uf T s ~dy 1 = ~dt> ' ~ J dz' 

Leading order corrections in x/l s will be retained only when multiplied by k y . That 
is, once the twisting nature of the perturbed structures arising due to the magnetic 
field shear has been accounted for by the frame transformation, only the leading effect 
of the background gradients, their drive of drift waves, is retained. The equation for 
h takes the form: 
dh 
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(75) 



The gyroaverages and Poisson bracket are now taken to represent their transformed 
values. The primes on the transformed variables are dropped from hereon. 

For clarity in identifying the effects of the flow shear, we consider the collisional 
limit of this system and so introduce the subsidiary ordering: 

v ^S> u, uj* , Vthk\\, Ufku, vk 2 p 2 . (76) 
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Here v is the ion-ion collision frequency, electron mass corrections are neglected, 
kp ~ O I yj u}jv\ <C 1 where the ion gyroradius is p = i>th/^ and oj* represents 
the drift frequencies associated with the background gradients, which arc identified 
in section 12.31 The distribution function is expanded in lo/v such that: h = 
hS ' + h( 1 ' + . . ., where the order of expansion is denoted by a superscript. We then 
proceed with the solution of the gyro-kinetic equation ([75]) order by order, taking all 
quantities on the right of ((76)) to be of the same order. 

Care must be taken when considering the ion self-collision operator. The Landau 
form, acting on a distribution function f (r,v,t), is [43] : 



<?(/./) = *£ 



dVll • 
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(77) 



where the relative velocity u = v — v' , the tensor U = (u 2 \ — uuj /u 3 , I is the identity 
tensor, Vi = e 4 In A/87regm 2 and In A is the Coulomb logarithm. The collision operator 
is Galilean invariant and for a Maxwellian, Fm (v) = (n /n 3 / 2 v 3 h ) exp (—v 2 /v 2 h ), 
C (Fm, Fm) = 0. Therefore the linearized collision operator appearing in (|75]) is: 

C l (h) = C(h,F M ) + C(F M ,h). (78) 

As discussed by Abel et al. [44], the velocity derivative appearing in the collision 
operator is taken at constant particle position, r, whilst the gyroaverage in (|75[) is 
taken at constant guiding centre position, R. As Fm does not vary on the scale of the 
gyroradius, this may be conveniently dealt with by using the Fourier representation 
of h, which is a function of guiding centre position, within the collision term: 
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Then, as kp <C 1: 
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Note that the gyroaverage eliminates terms of O (kp) in this expression. 
Now, to lowest order, ([75]) gives: 
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(81) 



This indicates [IS] that h^ 1 has the form of a perturbed Maxwellian in the guiding 
center variable R, with the perturbed density, Sn, temperature, ST and parallel 
velocity, 6Vf. 
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With all quantities evaluated at the guiding center position and the collisional 
contribution denoted as C^ 1 ' , the next terms in the expansion of (|75|l give: 
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The density, parallel velocity and energy moments of (|83l) give a set of equations 
describing the evolution of the parameters 5n, ST and SV\\ of the perturbed Maxwellian 
(|82[) . The collision term, C^\ has three contributions. The first, proportional to 
(k • p) C l (/i^ -*), is zero, given the form of h^ '. The second arises from terms of the 

form: C l (h k (k-p) ) and will not contribute to the moment equations due to the 

conservation properties of the collision operator: J (l, ii>|| , w 2 ) C l (g) dw = 0, where g 
is an arbitrary function. The remaining contribution is: 



Y j e lk - R (k- P C l {h k k-p)) R , 



(84) 



Defining: 



tp(w) = -h k (w)k- p = hk(w)w- cr, (85) 

b x k/il and using u • U = U • u = 0, the definition ((77)) shows that 



where cr = 


b 


x k/fl and 


reduces to: 




k 


with: 






[ d f ^ 


- 





w ■ crvj 



_d_ 

dw 



(scr ■ u 



dVU ■ 



ST k 



d i \ d i >\ 



dw 



dw' 



(86) 



n 



<T> 



(bu ■ cr + cru ■ b) 



25V\ 



\\k 



(87) 



■tii 



where we have used the sum and difference velocity variables: s = w + w' and 
u = w — w'. If we weight ([57]) by a function of velocity g (w\\,w±) and integrate 
over velocity space, integrating by parts noting that /m — > as \w\ — > 0, gives the 
resulting integral, M@, to be: 



fc = 


-5- /a/ a 


1(5- 


(scr 


N ST k 
u + au-s) 4 


(bu 



dg , dg 1 



dw 



dw' 



•U 



cru • b) 



26V, 



\\k 



: th 



l £) F »tii) 



where g' = g { w'„ , w' 



Now we take the density, parallel velocity and energy moments of (|83[) . which 
correspond to g = {l,wu, (y^w 2 — 3/2)}. For g = 1, the collisional contribution is 
seen to be zero by (|55|) . as required for particle conservation. To evaluate the collisional 
dissipation of the parallel velocity, note that the term in (|55|) proportional to ST is 
odd in s for g = w\\. So the dissipation is proportional to <$Vj|: 



'II 



^v th n 2 ^ 



e lk - u ki5V l{k 



= noi/Vi5V[|. (89) 

Similarly, the dissipation arising in the energy moment only retains a contribution 
from the term proportional to 5T in 



d d w 



w 

' th 



a (D 



AV2 Ujul y^ Jk-R k 2 ^± 
3 0F Vth « 2 Y ± T - 

3 n 



2T, 



XVI^T. 



(90) 
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Hence, the evolution equations of the perturbed variables are: 

d d\5n 3 1 f , 8n\ 1 d<f> efd d\, 

d-t +Uf dz)- Q + d-z 6V " + % T' - j = 1^8Jj + T l {d't + Uf d-z) ^ (91) 
d d \ „ r Id / m Sn ™\ 1 r , ,-,, i c 



8t •■ n 'Tz)^ + mTz[ T ^ + ^) + ^^-l±^ + V ^^ (92) 
d d \ ST 2 d r 1 ( ST] n d(j> ? ST 

m +U ^)ll + 3T Z SV ^^{^T-\ = ^B- /y + ^^ ^ 

where: 

<-*>-(!■§)& ^ 

For simplicity we now assume that the equilibrium ion and electron temperatures 
arc equal: Ti = T e = T. Upon imposing quasincutrality with the isothermal electron 
response given in (|67p. the electrostatic potential may be eliminated from these 
moment equations in favour of the perturbed density: 

Sn e __ e0 _ Sni _ /Sn _ ec/A. . . 

n T n \n T J ' 

This leaves a closed system for the evolution of the perturbed fluid variables: Sn, 
6V\\ and ST. We introduce the following normalizations, which are widely used in 
gyro-kinetic codes [38] : 

t = —i, (96) 

Cs 

■-. r=^-, *i= : -, (97) 

c s Ps ± Ps n p s 

where p s = c s /f2 is the sound Larmor radius. We also introduce an effective Mach 
number, associated with the convective effect of the perpendicular flow shear, defined 

as: 

M=^-. (98) 

Cs 

Dropping the tilde denoting the final coordinate transformation for convenience, the 
set of equations describing the evolution of a perturbation of the system finally takes 
the form: 



PsX, 


y = p s y, 


Z = IgZ, 


C s Ps' 


f= 5 Sk 

T Ps' 


Sn I 

n = 

n p 



dV _ 3 l s dn 
dz 8 l n dy ' 




(99) 
XJ-VlT, (101) 



3 l s dn 

8hdy 



where 



d 2 fd__ d_ 
dy 2 \ dx dy 



v i = 7TT+hr-^3- ( 102 ) 



Understanding the effect of sheared flow on microinstabilities 24 

and 

(y±,x±) : 



( Cs pi/i s y(c s pi/i s ), 

9 1\ [2l s n e 4 \nA 



40' 4 J V 3 87r 3 / 2 egT 2 



(103) 
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